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I review the field-theoretic renomalization group approach to quantum gravity, built 
around the existence of a non-trivial ultraviolet fixed point in four dimensions. I discuss 
the implications of such a fixed point, found in three largely unrelated non-perturbative 
approaches, and how it relates to the vacuum state of quantum gravity, and specifically 
to the running of G. One distinctive and well-known feature of the new fixed point is 
the emergence of a second genuinely non-perturbative scale, analogous to the scaling 
violation parameter in non-abelian gauge theories. I argue that it is natural to identify 
such a scale with the small observed cosmological constant, which in quantum gravity 
can arise as a non-perturbative vacuum condensate. I then show how the lattice cutoff 
theory of gravity can in principle provide quantitative predictions on the running of 
G, which can then be used to construct manifestly covariant effective field equations, 
and from there estimate the size of non-local quantum corrections to the standard GR 
framework. 
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1. Perturbative Non-renormalizability and Feynman Path Integral 

In a quantum theory of gravity the couphng constant is dimensionful, G ^ M^^"^, 
and within the standard perturbative treatment of radiative corrections one expects 
trouble in four dimensions, based on purely on dimensional grounds. The divergent 
one loop corrections are proportional to GA*^"^ where A is the ultraviolet cutoff, 
which then leads to a bad high momentum behavior, with an effective running 
Newton's constant 

G(/fc2)/G - l + ci(d) Gfc'^-2+ C»(G2) . (1) 

A more general argument for perturbative non-renormalizability starts by consid- 
ering the gravitational action with scalar curvature term R, which involves two 
derivatives of the metric. Then the graviton propagator in momentum space goes 
like and the vertex functions like . In c? dimensions each loop integral with 

involve a momentum integration d'^fc, so that the superficial degree of divergence T> 
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of a Feynman diagram with L loops is given by 

V^2 + {d-2)L , (2) 

independent of the number of external lines. One therefore concludes that for d > 2 
the degree of divergence for Einstein gravity increases rapidly with loop order L, and 
that the theory cannot be renormalized in naive perturbation theory. A consequence 
of the lack of perturbative renormalizability is the fact that new higher derivative 
counterterms arise to one-loop order [1] 

with even higher derivatives appearing at the next order. One concludes that the 
standard approach based on a perturbative expansion of the pure Einstein theory 
in four dimensions is not convergent; in fact it is badly divergent. 

A number of possible options have been proposed, the simplest of which is to 
just add the above higher derivative terms to the original action. The resulting ex- 
tended theory is perturbatively renormalizabe to all orders, but suffers potentially 
from unitarity problems. But these cannot be satisfactorily addressed in pertur- 
bation theory, as the theory is now asymptotically free in the higher derivative 
couplings and presumably exhibits a non-trivial vacuum. Alternatively, the gravity 
divergences can be cancelled by adding new unobserved massless particles and in- 
voking supersymmety; in fact it has been claimed recently that N = 8 supergravity 
might not be just renormalizable, but indeed finite to some relatively high loop 
order. The downside of this somewhat contrived approach is the artificial introduc- 
tion of a plethora of unobserved massless particles ofspinO, 1/2, 1,3/2, added to the 
original action in order to cancel the gravitational ultraviolet divergences. Finally, 
string theory solves the problem of ultraviolet divergences by postulating the exis- 
tence of fundamental extended objects, thus in part bypassing the requirement of 
supersymmetry and providing a natural cutoff for gravity, related to a fundamental 
string scale [2]. 

Nevertheless one important point that cannot be overlooked is the fact that 
in other field theories, which to some extent share with gravity the same set of 
ultraviolet problems (the non-linear sigma model is the most notable one, and the 
best studied case), the analogous result of Eq. ([T|) is in fact known to be incorrect. 
It appears as an artifact of naive perturbation theory, which in four dimensions 
does not converge, and seems to lead therefore to fundamentally incorrect answers. 
The correct answer in these models is found instead either by expanding around 
the dimension in which the theory is perturbatively renormalizable, or by solving it 
exactly in the large N limit and then computing 1/N corrections, or by solving it 
numerically on a lattice. The generic new feature in these models is the existence of 
a non-trivial fixed point of the renormalization group [3-7], which is inaccessible by 
perturbation theory in four dimensions, and radically alters the ultraviolet behavior 
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of the theory. The key question for gravity is therefore: are the uhraviolet problems 
just an artifact of a naive apphcation of perturbation theory in four dimensions, 
as clearly happens in other perturbatively non-renormalizable theories that also 
contain dimensionful couplings in four dimensions? 

In the following I will limit my discussion to the approach based on traditional 
quantum field theory methods and the renormalization group, applied to the Ein- 
stein action with a cosmological term, an avenue which in the end is intimately 
tied with the existence of a non-trivial ultraviolet fixed point in G in four dimen- 
sions. The nature of such a fixed point was first discussed in detail by K. Wilson 
for scalar and fermionic theories [3] , and the methods later applied to gravity in [6] , 
where they were referred to as asymptotic safety. As discussed above, it is fair to 
say that so far this is the only approach known to work in other not perturbatively 
renormalizable theories. 

If non-perturbative effects play an important role in quantum gravity, then one 
would expect the need for an improved formulation of the quantum theory, which 
does not rely exclusively on the framework of perturbation theory. After all, the 
fluctuating quantum metric field is dimensionless, and carries therefore no nat- 
ural scale. For the somewhat simpler cases of a scalar field and non-Abelian gauge 
theories a consistent non-perturbative formulation based on the Feynman path in- 
tegral has been used for some time, and is by now well developed. In a nutshell, 
the Feynman path integral formulation for quantum gravitation can be expressed 
by the functional integral formula 

J geometries 

Furthermore a bit of thought reveals that for gravity, to all orders in the weak field 
expansion, there is really no difference of substance between the Lorentzian (or 
pseudo-Riemannian) and the Euclidean (or Riemannian) formulation, which can 
be mapped into each other by analytic continuation. In the following therefore the 
Euclidean formulation will be assumed, unless stated otherwise. 

In function space one needs a metric before one can define a volume element. 
Therefore, following DeWitt, one first defines an invariant norm for metric defor- 
mations 

\\5gf = j d''x5g,,{x)G^^^^f[g{x)\59^p{x) , (5) 



" After QED, the second most accurate prediction of quantum field theory to date is for a per- 
turbatively non-renormalizable theory, the 0(N) non-linear cr-model in three dimensions, whose 
field theoretic treatment based on a non-trivial fixed point of the renormalization group, either on 
the lattice or in the continuum, eventually provides detailed predictions for scaling behavior and 
anomalous dimensions in the vicinity of the fixed point [5,8]. These have recently been verified ex- 
perimentally to high accuracy in a sophisticated space shuttle experiment [9] for critical superfuid 
Helium, whose order parameter corresponds to A'^ = 2 in the non-linear cr-model. 
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with the supermetric G given by the ultra-local expression 



G^^^-^f'[g{x)\ = \^)[g^^'^{x)g^l\x)+g'^P{x)g^'^{x)-\g^^^x)g'^f'{x)] (6) 

with A a real parameter, A ^ 2/d. The DeWitt supermetric then defines a suitable 
functional volume element a/G in four dimensions, 

/ [dg^A = / n n '^9,.{x) . (7) 

The Euclidean Feynman path integral for pure Einstein gravity with a cosmological 
constant term is then written as 



j [dgij^i,] exp |-Ao J dx^ + J dx^R^ . 



(8) 



An important aspect of this path integral is connected with the global scaling prop- 
erties of the action and the measure [10]. First one notices that in pure Einstein 
gravity with a bare cosmological constant term 

^ = -16^^'' + '°^ 
one can rescale the metric by g^i, = lv with uj a constant, giving 

C = ~--^oj''^'-'^R' + XoLo^/'^ . (10) 
IdttGo 

This can then be interpreted as a rescaling of the two bare couplings Go — >■ 
^-d/2+i(^^^ Ao — )• Aqw'^/^, leaving the dimensionless combination GqAq"^ un- 
changed. Therefore only the latter combination has physical meaning in pure gravity. 
In particular, one can always choose the scale a; = Ag ^^'^ , so as to adjust the vol- 
ume term to have a unit coefhcient. The implication of this last result is that pure 
gravity only contains one bare coupling Go, besides the ultraviolet cutoff A needed 
to regulate the quantum theory. 



2. Gravity in 2 + e Dimensions and Non- Trivial UV Fixed Point 

In two dimensions the gravitational coupling becomes dimensionless, G ~ A^"'', 
and the theory appears perturbatively renormalizable. In spite of the fact that the 
gravitational action reduces to a topological invariant, it is meaningful to attempt 
to construct, in analogy to what was suggested originally by Wilson for scalar field 
theories, the theory perturbatively as a double series in e = d — 2 and G. The 2 + e 
expansion for pure gravity then proceeds as follows [6,11]. First the gravitational 
part of the action 

C = -^^^R + Xo^, (11) 

with G now dimensionless and fj, an arbitrary momentum scale, is expanded by 
setting 

9ixu g^v = g^v + h^v , (12) 



February 3, 2010 18:11 



WSPC - Proceedings Trim Size: 9.75in x 6.5in mgl2 



5 

where g^i, is the classical background field, and ft.^^ a small quantum fluctuation. 
The quantity C in Eq. (|lll) is naturally identified with the bare Lagrangian, and 
the scale fi with a microscopic ultraviolet cutoff A, corresponding to the inverse 
lattice spacing in the lattice formulation. After the quantum fluctuations in /i^^ are 
integrated out and the cosmological constant term gets rescaled, one obtains the 
following result for the renormalization group beta function for G: with Ns scalar 
fields and Np Majorana fermion fields the result to two loops reads [12] 

n^G = 13(G) =eG - PoG^ ~ (3iG^ + OiG'^, G^e, G^e^) , (13) 

with /3o = I (25 - c) and /3i = ^ (25 - c), and c = A^s + A^f/2. The physics of this 
result is contained in the fact that the gravitational /3-function determines the scale 
dependence of Newton's constant G, and has the shape shown in Fig. 1. 



13(G) 




Figure 1. Renormalization group /3-function for gravity close to two dimensions. The arrows 
indicate the coupling constant flow towards increasingly larger distance scales. 

A closer examination of the result shows that depending on whether one is on 
the right (G > Gc) or on the left (G < Gc) of the non-trivial ultraviolet fixed point 
at 

d~2 



Gr 



/So 



o{{d-2y) 



(14) 



(with Gc positive, provided one has c < 25) the coupling will either flow to in- 
creasingly larger values of G, or flow towards the Gaussian flxed point at G = 0, 
respectively. Furthermore the running of G as a function of the sliding momentum 
scale n = k can be obtained by integrating Eq. p^ . and one has to lowest order 

(d-2)/2 



G(fc^) 



Gc 



, m 
1 ± ^0 ( ^ 



(15) 



with Co a positive constant, and m a new nonperturbative scale. As in non-abelian 
gauge theories and QCD, this last quantity arises naturally as an integration con- 
stant of the renormalization group equations. The choice of + or — sign is then 
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determined from whether one is to the left (+), or to right (-) of Gc, in which 
case the effective G{k'^) decreases or, respectively, increases as one flows away from 
the ultraviolet fixed point towards lower momenta, or larger distances. Physically 
therefore the two solutions represent a gravitational screening (G < Gc), and a 
gravitational anti-screening (G > Gc) situation [10]. Finally, at energies sufficiently 
high to become comparable to the ultraviolet cutoff, the gravitational coupling G 
eventually flows towards the ultraviolet fixed point G{k^) ^^^^f^i G(A), where 
G(A) is the coupling at the cutoff scale A, to be identified with the bare or lattice 
coupling. 

One message therefore is that the quantum corrections involves a new physical, 
renormalization group invariant scale ^ = 1/m-, which cannot be fixed perturbatively 
and whose size determines the scale for the new quantum effects. In terms of the 
bare couphng G(A), it is given by 

/ /■G(A) dG' \ 
r'-m = A™.Aexpf-y j . (16) 

The constant Am on the r.h.s. of Eq. (fTB)) cannot be determined by perturbation 
theory; it needs to be computed by non-perturbative lattice methods. 

At the fixed point G = Gc the theory is scale invariant by definition; in statis- 
tical field theory language the fixed point corresponds to a phase transition. In the 
vicinity of the fixed point one can write 

/3(G) ^ /3'(Gc) (G - Gc) + 0((G - Gc)^) . (17) 

If one defines the exponent v by fi'{Gc) — — 1/i^, then from Eq. ([TE)) one has by 
integration 

m - A- An|G(A)-Gcr , (18) 

G — c 

with u the correlation length exponent. To two loops the results of [12] imply 

15 



25 - c 



.£2 



(19) 



which gives, for pure gravity without matter (c = 0) in four dimensions, to lowest 
order the scaling exponent ~ 2, and w 4.4 at the next order. The key 
question raised by these 2 + e perturbative calculations is therefore: what remains 
of the above phase transition in four dimensions, how are the two phases of gravity 
characterized non-perturbatively, and what is the value of the exponent v deter- 
mining the running of G in the vicinity of the fixed point in four dimensions ? To 
answer this question in a controlled way would seem to require the introduction of a 
non-perturbative regulator, based on the lattice formulation (since no other reliable 
non-perturbative regulator for field theories is known to date). 
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3. Lattice Regularized Quantum Gravity 

On the lattice the infinite number of degrees of freedom in the continuum is re- 
stricted, by considering Riemannian spaces described by only a finite number of 
variables, to the geodesic distances between neighboring points. Such spaces are 
taken to be flat almost everywhere, and referred to as piecewise linear. The elemen- 
tary building blocks for d-dimensional space-time are then simplices of dimension d. 
A 0-simplex is a point, a 1-simplex is an edge, a 2-simplex is a triangle, a 3-simplex 
is a tetrahedron. A d-simplex is a d-dimensional object with d + 1 vertices and 
d{d + l)/2 edges connecting them [13]. 

The geometry of the interior of a d-simplex is assumed to be flat, and is therefore 
completely specifled by the lengths of its d{d + l)/2 edges. When focusing on one 
such n-simplex, it is convenient to label the vertices by 0, 1, 2, 3, ... , n and denote 
the square edge lengths by Iqi = I^q, ... , /q„. The simplex can then be spanned 
by the set of n vectors ei, ... e„ connecting the vertex to the other vertices. To 
the remaining edges within the simplex one then assigns vectors Cij = Ci — ej with 
1 ^ * < j ^ ^- Within each n-simplex one can define a metric gij{s) = • ej, and 
then in terms of the edge lengths lij — \ei — ej\ the metric is given by 

= h (^o. + 'o, -^',) • (20) 

The volume of a general n-simplex can be found by the n-dimensional generalization 
of the well-known formula for a tetrahedron, namely 

F„(s) = l^det<?,,(s) . (21) 

In a piecewise linear space curvature is detected by going around elementary loops 
which are dual to a (d — 2)-dimensional subspace. From the dihedral angles 9{s, h) 
associated with the faces of the simplices meeting at a given hinge h one computes 
the deficit angle 5{h), defined as [13] 

5{h) ^ 27T-J2 (^is,h) , (22) 

sDh 

where the sum extends over all simplices s meeting on h. It then follows that the 
deficit angle ^ is a measure of the local curvature at h. By considering rotation 
matrices around a hinge one can obtain an expression for the Riemann tensor at 
the hinge h 

R^.xaih) = -r^C^p.(/i)C^A.(/^) (23) 

which is expected to be valid in the limit of small curvatures, with Ac{h) the area of 
the loop entangling the hinge h, and U a bivector describing the hinge's orientation. 
From the expression for the Riemann tensor at a hinge given in Eq. (|23p one obtains 
by contraction the Ricci scalar R(h) — 2 5{h) / Ac{h), and the continuum expression 
^/g R is then obtained by multiplication with the volume element V{h) associated 
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with a hinge. The curvature and cosmological constant terms then lead to the 
combined Regge lattice action 

^ Xo J2 V('^-k Y: S,vt'-'\ (24) 

simpliccs s hinges h 

One key aspect of this formulation is the local gauge invariance of the lattice action, 
in analogy to the local gauge invariance of the Wilson action for gauge theories. Al- 
ready on a flat 2-d lattice it is clear that one can move around a point on a surface, 
keeping all the neighbors fixed, without violating the triangle inequalities, and leav- 
ing local curvature invariants unchanged. In d dimensions this transformation has d 
parameters and is an exact invariance of the action. When space is slightly curved, 
piecewise linear diffeomorphisms can still be defined as the set of local motions of 
points that leave the local contribution to the action, the measure and the lattice 
analogs of continuum curvature invariants unchanged. In the limit when the number 
of edges becomes very large one expects the full continuum diffeomorphism group 
to be recovered [14]. 

In order to write down a lattice path integral, one needs, besides the action, a 
functional measure. As the edge lengths hj play the role of the continuum metric 
gnv{x), one expects the discrete measure to involve an integration over the squared 
edge lengths. After choosing coordinates along the edges emanating from a vertex, 
the relation between metric perturbations and squared edge length variations for a 
given simplex based at in d dimensions is from Eq. pop 



59.,,{f) = ^i6ll, + Sll-5lf^) . (25) 

For one d-dimensional simplex labeled by s the integration over the metric is thus 
equivalent to an integration over the edge lengths, and one has the identity 

. d(d_i) ci(d+l)/2 

-y^detg,,(s)) n dg^jis) = (-i)^ [V^l^)] J] '^^k- (26) 

^ i>j k=l 

Indeed there are d{d + l)/2 edges for each simplex, just as there are d{d + l)/2 
independent components for the metric tensor in d dimensions. In addition, a certain 
set of simplicial inequalities need to be imposed on the edge lengths. These represent 
conditions conditions on the edge lengths kj such that the sites i can be considered 
as vertices of a d-simplex embedded in flat d-dimensional Euclidean space. After 
summing over all simpliccs one derives what is regarded as the lattice functional 
measure representing the continuum DeWitt measure in four dimensions 

J Jo 

Here & is a (step) function of the edge lengths, with the property that it is equal 
to one whenever the triangle inequalities and their higher dimensional analogs are 
satisfied and zero otherwise. The lattice action of Eq. ([24|) for pure four-dimensional 
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Euclidean gravity then leads to the regularized lattice functional integral [7] 

Ziatt = jW] expj-Ao^K + fc^(5,,A,,| , (28) 

where, as customary, the lattice ultraviolet cutoff is set equal to one (i.e. all length 
scales are measured in units of the lattice cutoff). Furthermore, Ao sets the overall 
scale, and can therefore be set equal to one without any loss of generality, according 
to the scaling arguments presented before. The lattice partition function Ziatt should 
then be compared to the continuum Euclidean Feynman path integral for pure 
gravity of Eq. ©. 

In closing we note that what makes the Regge theory stand out compared to 
other possible discretization of gravity is the fact that it is the only lattice theory 
known to have the correct spectrum of continuous excitations in the weak field limit, 
i.e. transverse traceless modes, or, equivalently, helicity-two massless gravitons. In- 
deed one of the simplest possible problems that can be treated in lattice quantum 
gravity is the analysis of small fluctuations about a fixed flat simplicial background. 
In this case one finds that the lattice graviton propagator in a De Donder-like gauge 
is identical to the continuum expression [7]. 



4. Strongly Coupled Gravity and Gravitational Wilson Loop 

As in non-abelian gauge theories, important information about the non-perturbative 
ground state of the theory can be gained by considering the strong coupling limit. 
In lattice gravity an expansion can be performed for large G or small k = I/SttG, 
and the resulting series is in general expected to be useful up to some k — kc, where 
kc is the lattice critical point, at which point the partition function Z eventually 
develops a singularity. One starts from the lattice regularized path integral with 
action Eq. ([24|) and measure Eq. ([27|) . The four-dimensional Euclidean lattice action 
usually contains a cosmological constant and scalar curvature term as in Eq. (p4|) . 

Ilatt ^ \Y,Ml^)- kJ2Sk{P)Ml^) ■ (29) 

h h 

The action only couples edges which belong either to the same simplex or to a set of 
neighboring simplices, and can therefore be considered local just like the continuum 
action; it leads to the lattice partition function defined in Eq. (^5]) . When doing 
an expansion in the kinetic term proportional to fc, it is convenient to include the 
A-term in the measure, and Ziatt can then be expanded in powers of k, 

Ziattik) = / dfi{l') e''^-'-^- = E V" / '^^(^'^ (T.^'^^^] ■ (30) 

n=0 '"'■J \ h ) 

Z(k) — 'Y^=Q '^n k"' is analytic at A: = 0, so this expansion is well defined up to the 
nearest singularity in the complex k plane. 
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In the gravity case the analogs of the gauge variables of Yang-Mills theories 
are given by the connections, so it is natural when computing the gravitational 
Wilson loop [15] to look for a first order formulation of Regge gravity [16]. For each 
neighboring pair of simplices s, s + 1 one can associate a Lorentz transformation 
R^^,{s,s + 1), and one then might want to consider a near-planar closed loop C, 
such as the one shown schematically in Fig. 2. Along a closed loop the overall rotation 
matrix is given by 

1 



s CC 



(31) 




Figure 2. Gravitational analog of the Wilson loop. A vector is parallel-transported along the 
larger outer loop. The enclosed minimal surface is tiled with parallel transport polygons, 
here chosen to be triangles for illustrative purposes. For each link of the dual lattice, the 
elementary parallel transport matrices R(s,s') are represented by arrows. 

In the quantum theory one is interested in averages of the above product of 
rotations along a given path. If the curvature of the manifold is small, then clas- 
sically the expression for the rotation matrix R(C) associated with a near-planar 
loop can be re-written in terms of a surface integral of the large-scale Riemann 
tensor, projected along the surface area element bivector A"'^(C) associated with 
the orientation of the loop, 

i / R-.^pA^P{C) , 

Rf'.iC) ^ e Js ^ . (32) 

L J 

Thus a direct calculation of the quantum Wilson loop could in principle provide a 

way of determining the effective curvature on very large distance scales, even in the 

case where short distance fluctuations in the metric may be significant. 

A detailed lattice calculation in the strong coupling limit then gives the following 

result. First one defines the lattice Wilson loop as 

W{C) = < Tr[{Uc + e h) Ri R2 Rn] > ■ (33) 

where the Rfs are the rotation matrices along the path and the factor {Uc + £^4) 
contains some average direction bivector Uc for the loop, which is assumed to be 
close to planar. Then for sufficiently strong coupling one can show that one naturally 
obtains an area law, which here we express as 

W{C) ^ exp{-Ac/e) (34) 
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where ^ is the gravitational correlation length. The appearance of f follows from 
dimensional arguments, given that the correlation length is the only relevant length 
scale in the vicinity of the fixed point; the results can thus be considered analogous 
to the well-known behavior for the Wilson loop in non-abelian gauge theories [17]. In 
the actual calculation the rapid decay of the quantum gravitational Wilson loop as 
a function of the area can be seen as a general consequence of the assumed disorder 
in the uncorrelated fluctuations of the parallel transport matrices R(s, s') at large 
G. A careful identification of (a suitable trace of) Eq. (|32]) with the expression in 
Eq. (j34p . and in particular the comparison of the two area-dependent terms, then 
yields the following estimate for the macroscopic, large scale, average curvature in 
the large G limit 

R - i/e , (35) 

where ^ is the quantity in Ec^. ()16|) . An equivalent way of phrasing the last result 
is the suggestion that l/^'^, where ^ is the renormalization group invariant grav- 
itational correlation length of Eq. (fT6|) . should be identified, up to a constant of 
proportionality of order one, with the observationally determined, large scale cos- 
mological constant A. 

5. Nonperturbative Gravity 

The exact evaluation of the lattice functional integral for quantum gravity by nu- 
merical methods allows one, in principle, to investigate a regime which is generally 
inaccessible by perturbation theory: where the coupling G is strong and quantum 
fluctuations in the metric are large. The hope is, in the end, to make contact with 
the analytic results obtained in the 2 + e expansion, and determine which scenar- 
ios are physically realized in the lattice regularized model. The main question one 
would therefore like to answer is whether there is any indication that the non-trivial 
ultraviolet fixed point scenario is realized in the lattice theory, in four dimensions. 
This would imply, as in the non-linear sigma model and similar models, the exis- 
tence of at least two physically distinct phases, and associated non-trivial scaling 
dimensions. A clear physical characterization of the two gravitational phases would 
also allow one, at least in principle, to decide which phase, if any, could be realized 
in nature. As discussed below, the lattice continuum limit is taken in the vicinity 
of the fixed point, so close to it is perhaps the physically most relevant regime. 

At the next level one would hope to be able to establish a quantitative connection 
with the continuum perturbative results, such as the 2 + e expansion discussed 
earlier. Since the lattice cutoff and the method of dimensional regularization cut the 
theory off in the ultraviolet in rather different ways, one needs to compare universal 
quantities which are cutoff-independent. An example is the critical exponent v, as 
well as any other non-trivial scaling dimension that might arise. One should note 
that within the 2 + e expansion only one such exponent appears, to all orders in 
the loop expansion, as i'^^ = —l3'{Gc)- Therefore one central issue in the four- 
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dimensional lattice regularized theory is the value of the universal scaling exponent 
V [10,18] [in Eqs. ([T5 |) . ([T^ and ] . 

The starting point is again the lattice regularized path integral with action as 
in Eq. (j24[) and measure as in Eq. (j27p . Among the simplest quantum mechanical 
averages that one can compute is one associated with the local curvature, 

But the curvature associated with this quantity is one that would be detected when 
parallel-transporting vectors around very small infinitesimal loops. Furthermore 
when computing correlations in quantum gravity new subtle issues arise, due to 
the fact that the physical distance between any two points x and y 

dix, y\g) = min dr ^/^^^(O^f^ (37) 

? Jt(x) 

is a fluctuating function of the background metric g^j_y{x). Consequently physical 
correlations have to be defined at fixed geodesic distance d, as in the following 
connected correlation between observables O 

<Jdx j dy^0(x)^0{y)8{\x-y\-d)>, . (38) 

Based on general arguments one expects such correlations to either follow a power 
law decay at short distances, or an exponential decay characterized by a correlation 
length ^ at larger distances 

<^gO{x)^0{y)5{\x-y\-d)>,^ e-''/^ . (39) 

In practice such correlations at fixed geodesic distance are difhcult to compute 
numerically, and therefore not the best route to study the critical properties of 
the theory. But scaling arguments allow one to determine the scaling behavior of 
correlation functions from critical exponents characterizing the singular behavior of 
the free energy F{k) = —{1/V) In Z and of various local averages in the vicinity of 
the critical point. In general a divergence of the correlation length ^ 

m ^-^^ \kc-kr (40) 

signals the presence of a phase transition, and leads to the appearance of a non- 
analyticty in the free energy F{k). One way to determine v is from the curvature 
fluctuation, for which one can show 

Xn{k) ^^^^ A^^ \k,-k\'(^-'''^) . (41) 

From such averages and fluctutations one can therefore, in principle, extract the 
correlation length exponent v of Eq. ()40p , without having to compute an invariant 
correlation function at fixed geodesic distance. 

In general for the measure in Eq. (|27p one finds a well behaved ground state only 
for k < kc [10]. The system then resides in the 'smooth' phase, with an effective 
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dimensionality close to four. On the other hand, for k > kc the curvature becomes 
very large and the lattice collapses locally into degenerate configurations with very 
long, elongated simplices. This last phenomenon is usually interpreted as a lattice 
remnant of the conformal mode instability of Euclidean gravity. 

There are a number of ways by which the critical exponents can be determined 
to some accuracy from numerical simulations, and it is beyond the scope of this short 
review to go into more details; one obtains eventually in = 4 /cc — 0.0636 and 
I' ~ 0.335, which suggests 

1/^1/3 (42) 

for pure quantum gravity in four dimensions [18]. Note that at the critical point the 
gravitational coupling is not weak, since Gc ~ 0.626 in units of the ultraviolet cutoff. 
Fig. 3 shows a comparison of the critical exponent v obtained by three independent 
methods, namely the original lattice result in d = 2,3,4 [18], the recent 2 + e 
expansion (to one and two loops) [12], and the even more recent renormalization 
group truncation method [19-21]. 
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Figure 3. Universal renormalization group scaling exponent l/z^ of Eq. (|17|) . computed 
in the lattice theory of gravity (large dots for d — 2, d = 3 and d = 4, and continous 
interpolating line) [18] . For comparison the 2 + e result is shown (two lower curves, to one 
(lower) and two (upper) loops) [12], as well as the recent Einstein-Hilbert truncation result 
(smaller dots and connecting dashed line) [21]. The abscissa is, for convenience, a variable 
related to the space-time dimension d through z = (d — 2)/(d — 1), which maps d = oo 
(where it is known that u = 0) to z = 1. Note that for a scalar theory one has 1/v = 2 for 
d > 4 [3]. 
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6. Renormalization Group, Lattice Continuum Limit and the 
Running of G 

The lattice theory points to the existence of a phase transition in pure quantum 
gravity, with a divergent correlation length in the vicinity of the critical point, as in 
Eq. ()40[) , which can be re- written in terms of the inverse correlation length to = 1/^ 

= TO = AA„,\kc - fcr . (43) 

In the above expression the correct dimensions have been restored, by inserting 
explicitly on the r.h.s. the ultraviolet cutoff A. Here k and fee are dimensionless 
quantities, corresponding to bare microscopic couplings at the cutoff scale, k = 
k{A) = l/87rG(A). is related to in Eq. (gD]) by Am = A~^. It is worth 
pointing out that the above expression for m{k) is identical in structure to the 2 + e 
result for continuum gravity, Eq. (fT8)) . 

Then the lattice continuum limit corresponds to a large cutoff limit, taken at 
fixed TO or ^, 

A — >■ 00 , k ^ kc , TO fixed , (44) 

which shows that the continuum limit is in fact reached in the vicinity of the ul- 
traviolet fixed point, k ^ kc- In practice, since the cutoff ultimately determines 
the physical value of Newton's constant G, the cutoff A cannot be taken to oo, 
and it persists as a fundamental scale in the theory. A very large value will suffice, 
A~^ ^ 10~^'^CTO, for which it will still be true that ^ ^ A~^, which is all that is 
required for the continuum limit. 

In order to discuss the renormalization group behavior of the coupling in the 
lattice theory it is convenient to re-write the result of Eq. (j43p directly in terms of 
Newton's constant G as 

with the dimensionless constant cq related to Am by Am — l/{cokc)'^. The above 
expression only involves the dimensionless ratio G{A)/Gc, which is the only relevant 
quantity here. From the knowledge of the dimensionless constant Am in Eq. (j43p one 
can estimate from first principles the value of cq in Eq. ([45|) and later in Eq. (j49|) . 
Lattice results for the correlation functions at fixed geodesic distance give a value 
for Am ~ 0.72 with a significant uncertainty, which, when combined with the values 
fcj, ~ 0.0636 and v ~ 0.335 given above, gives cg — l/ik^All") ^ 42. Then the 
renormalization group invariance of to = requires that the running gravitational 
coupling G(/i) varies in the vicinity of the fixed point in accordance with the above 
equation, with A — >■ ^, where is an arbitrary momentum scale, 
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The latter is equivalent to the renormalization group invariance requirement 

^^rniiiMp)) =Q (47) 

provided G{^) is varied in a specific way. Thus Eq. ()47|) can be used to obtain a 
Callan-Symanzik /3-function for the coupling G{^) in units of the ultraviolet cutoff, 

/i^G(A*) = , (48) 

with 13(G) given in the vicinity of the non-trivial fixed point, from Eq. (j46p . by 
/3(G) ^G-^Gc Gc)- Or one can obtain the scale dependence of the gravita- 

tional coupling directly from Eq. (j46p . which then gives 

G(Ai) = G, [l + co(mV/i')'/"^ + 0((mVfc2)^)] (49) 

in the physical anti-screening phase. Again, this last expression can be compared 
directly to the lowest order 2 -I- e result of Eq. (ITSt . The physical dimensions of G can 
be restored, by multiplying the above expression on both sides by the ultraviolet 
cutoff A, if one so desires. Physically the above lattice result implies anti-screening: 
the gravitational coupling G increases slowly with distance. 



7. Curvature Scales and Gravitational Condensate 

The renormalization group running of G(/i) in Eq. (|49)) involves an invariant scale 
^ = 1/m. At first it would seem that such a scale could take any value, including 
a very small one based on the naive estimate ^ ^ Ip - which would then preclude 
any observable quantum effects in the foreseeable future. But the results from the 
gravitational Wilson loop at strong coupling would suggest otherwise, namely that 
the non-perturbative scale ^ is in fact related to macroscopic curvature. From astro- 
physical observation the average curvature is very small [22] , so one would conclude 
that ^ has to be very large and possibly macroscopic, 

Aobs - ^ (50) 

with Xobs the observed small but non- vanishing scaled cosmological constant. A 
further indication that the identification of the observed cosmological constant with 
a mass-like - and therefore renormalization group invariant - term might make sense 
beyond the weak field limit can be seen, for example, by comparing the structure 
of the three classical field equations 

Rfiu — \ gfj.u R + X g^j.v — SttG T^i, 

for gravity, QED (made massive via the Higgs mechanism) and a self-interacting 
scalar field, respectively. Nevertheless it seems so far that the strongest argument 
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suggesting the identification of the scale ^ with A is derived from the calculation of 
the gravitational Wilson loop at strong [16]. 

This relationship, taken at face value, implies a very large, cosmological value 
for ^ ^ lO^^CTO, given the present bounds on Xphys- Thus a set of modified Einstein 
equations, incorporating the quantum running of G, would read 

Rfiv - R + '>^gij,u ^ 8TrG{p)T^^ (52) 

with A ~ and G(/i) on the r.h.s. scale-dependent, in accordance with Eq. (j49|) . 
The precise meaning of C?(/Lt) in a covariant framework is given below. 



8. Effective Covariant Field Equations 



The result of Eq. (|49p implies a running gravitational coupling in the vicinity of 
the ultraviolet fixed point, with m = 1/^, cq > and v ~ 1/3. Since f is expected 
to be very large, the quantity Gc in the above expression should now be identified 
with the laboratory scale value -\/G7 ~ 1.6 x 10~^^cto. The effective interaction in 
real space is then obtained by Fourier transform, but since the above expression 
is singular as — > 0, the infrared divergence needs to be regulated, which can 
be achieved by utilizing as the lower limit of momentum integration m = 1/^. A 
properly infrared regulated version of the above would read 



Gc 



1 



Co 



fc2 



+ 



(53) 



Then at very large distances r 3> the gravitational coupling is expected to ap- 
proach the finite value Goo = {I + cq + ■ ■ ■) Gc- 

The first step in analyzing the consequences of a running of G is to re-write 
the expression for G{k^) in a coordinate-independent way, for example by the use 
of a non-local Vilkovisky-type effective action. Since in going from momentum to 
position space one usually employs — > — □, to obtain a quantum-mechanical 
running of the gravitational coupling one has to make the replacement G — > G{0). 
Therefore from Eq. ((49)) one obtains 

l/2i/ 

1 + Co ' 



G(n) = Gc 



1 



(54) 



and the running of G is expected to lead to a non-local gravitational action, for 
example of the form 

l/2iy 

1 - Co 



I 



16nG 



1 



R 



(55) 



Due to the appearance of a fractional exponent, the covariant operator appearing 
in the above expression has to be suitably defined by analytic continuation. The 
latter can be done, for example, by computing for positive integer n and then 
analytically continuing to rt — > —1/21^. 
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Alternatively, had one not considered the action of Eq. (1551) as a starting point 
for constructing the effective theory, one would naturally be led (as suggested by 
Eq. ([M))) to consider instead the following effective field equation 

-Rm. - ^ff^i^i? + Ag^, = 87rG(n)r^, , (56) 

the argument again being the replacement G — > G(n) in the classical Einstein field 
equations. Being manifestly covariant, these expressions at least satisfy some of the 
requirements for a set of consistent field equations incorporating the running of G. 

The effective field equations of Eq. ([56l) can in fact be re-cast in a form very 
similar to the classical field equations but with a T^i, = [G(n)/G] T^i, defined as an 
effective, or gravitationally dressed, energy-momentum tensor. Just like the ordinary 
Einstein gravity case, in general T^^ might not be covariantly conserved a priori, 
V'^ Tfj^u 7^ 0, but ultimately the consistency of the effective field equations demands 
that it be exactly conserved, in consideration of the Bianchi identity satisfied by 
the Riemann tensor. In this picture therefore the running of G can be viewed as 
contributing to a sort of " vacuum fluid" , introduced in order to to account for the 
gravitational vacuum polarization contribution. 



9. Static Isotropic Solutions 

One can show that the quantum correction due to the running of G can be described 

- at least in the non-relativistic limit of Eq. ([53)) when applied to Poisson's equation 

- in terms of a vacuum energy density Pmif), distributed around the static source 
of strength M in accordance with 

Pr,i{r) = -^c^coMm^ {mr)-i^^-i'> Ki(^r^_3_^{mr) (57) 

and with C;. = 25(5-^)/0Fr(^), and 

4tt r^drp,n{r) = cq M . (58) 
Jo 

More generally in the fully relativistic case, after solving the covariant effective field 
equations with G(n) for i/ = 1/3 one finds in Schwarzschild coordinates, and in the 
limit r > 2MG, 

A-\r) == B(r) = 1 + — ^— In mr + ... (59) 

r 3 TT 

The last expressions for A{r) and B{r) are therefore consistent with a gradual slow 
increase in G with distance, in accordance with the formula 

G ^ G(r) ^ G (1 + rMn + . . . ) (60) 

\ 3 vr TO^ ) 

in the regime r ^ 2M G. The last result is in some ways reminiscent of the QED 
small-r result 

Q ^ Q{t) ^ Q (l + + ..\ (61) 
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In the gravity case, the correction vanishes as r goes to zero: in this hmit one is 
probing the bare mass, unobstructed by its virtual graviton cloud. In some ways the 
running G term acts as a local cosmological constant term, for which the r depen- 
dence of the vacuum solution for small r is fixed by the nature of the Schwarzschild 
solution with a cosmological constant term. One could in fact wonder what these 
solutions might look like in d dimensions, and after some straightforward calcula- 
tions one finds that in d > 4 space-time dimensions a solution to the effective field 
equations can only be found if in Eq. (|54p v ~ l/(<i — 1) exactly [23]. 

10. Cosmological Solutions 

A scale dependent Newton's constant is expected to lead to small modifications 
of the standard cosmological solutions to the Einstein field equations. Here I will 
summarize what modifications are expected from the effective field equations on 
the basis of G(n), as given in Eq. ((54|). which itself originates in Eqs. ((53|) and (j49|) . 
The starting point are the quantum effective field equations of Eq. ([56| . with G'(n) 
defined in Eq. ([54| . In the Friedmann- Robertson- Walker (FRW) framework these 
are applied to the standard homogeneous isotropic metric 

dT^ = df - a" [t) I + {d9^ + sin^ d^^) | ^q^) 

It should be noted that there are in fact two related quantum contributions to the 
effective covariant field equations. The first one arises because of the presence of 
a non- vanishing cosmological constant A ~ caused by the non-perturbative 

vacuum condensate of Eq. ([50]). As in the case of standard FRW cosmology, this 
is expected to be the dominant contributions at large times t, and gives an expo- 
nential (for A > 0), or cyclic (for A < 0) expansion of the scale factor. The second 
contribution arises because of the explicit running of G{0) in the effective field 
equations. The next step is therefore a systematic examination of the nature of 
the solutions to the full effective field equations, with G{0) involving the relevant 
covariant d'Alembertian operator 

□ = 5^'^V,,V, (63) 
acting on second rank tensors as in the case of T^i, , 

(V^Tap) = d^^I^aP ~ ^If^IxaP ^ ^a^IvXP ^ ^^Pf^haX ■ (64) 

To start the process, one assumes for example that T^ti/ has a perfect fluid form, 
for which one obtains the action of on T^^^ and then analytically continues 
to negative fractional values of n = —1/2;^. Even in the simplest case, with G{U) 
acting on a scalar such as the trace of the energy-momentum tensor T\ , one finds 
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for the choice p{t) — p^t^ and a{t) = VQt" the rather unwieldy expression 



hp(i)] ^ 4" (-1) 



n+l 



r(f 



p ^ /3+3a+l 



1 



f3-2n 



(65) 



with integer n then analyticahy continued to n — — j^, with = 1/3. 

A more general calculation shows that a non- vanishing pressure contribution is 
generated in the effective field equations, even if one initially assumes a pressureless 
fluid, p{t) — 0. After a somewhat lengthy computation one obtains for a universe 
filled with non-relativistic matter (p—O) the following set of effective Friedmann 
equations, namely 



k 



3 

BttG r 



1 



1 

ct it/0'/- 



Pit) 



(66) 



for the tt field equation, and 

k a^it) 2a(i) 



SttG 



Ct it/0 



+ 



Pit) + A (67) 



for the rr field equation. In the above equations the running of G appropriate for 
the Robertson- Walker metric is given by 



Git) = G 



1 



Ct 



+ 



(68) 



with Ct of the same order as cq of Eq. . Note that it is the running of G that 
induces an effective pressure term in the second (rr) equation, corresponding to 
the presence of a relativistic fluid arising from the vacuum polarization contribu- 
tion. The second important feature of the new equations is an additional power-law 
acceleration contribution, in addition to the standard one due to A. 



11. Quantum Gravity and Cosmological Density Perturbations 

Besides the cosmic scale factor evolution and the static isotropic solutions just dis- 
cussed, the running of G(n) also affects the nature of matter density perturbations 
on very large scales. In discussing these effects, it is customary to introduce a per- 
turbed metric of the form 

dT^ = d£- ~ (% -I- /ly ) dx'^dx^ (69) 

with aii) the unperturbed scale factor and hijix,t) a small metric perturbation. 
The next step is to determine the effects of the running of G on the relevant matter 
and metric perturbations, again by the use of the modified field equations. For 
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sufficiently small perturbations, one can expand G{0) appearing in the effective 
covariant field equations in powers of the metric perturbation hij as 



G{a) = Go 



(70) 



It is also customary to expand the density, pressure and metric trace perturbation 
modes in spatial Fourier components 

Sp{x,t)^Sp{t)e'''-^ 5p{x,t) = 5p{t)e'^-^ h{x,t) = h{t) e'^-^ . (71) 

Normally the Einstein field equations R^v — ^ff/if dixu = SttG T^h are given to 
first order in the small perturbations by 

^^h{t) = %T,Gp{t)5(t) 

h{t) + 3 ^ /i(t) = -247r G w p{t)8{t) (72) 

a\t) 

with 5{t) = Sp{t)/p{t) and w = for non-relativistic matter, yielding then a single 
equation for the trace of the metric perturbation h{t), 

h{t) + 2 '^h{t) = -87rG(l + 3 w)p{t)5{t) . (73) 

Combined with the first order energy conservation {1 + w) h{t) = 5{t), this 
then gives a single equation for the density contrast S{t), 

6{t) + 2 ^ 5{t) -AttG p{t) 5{t) = . (74) 

In the case of a running G(n) these equations need to be re-derived from the effective 
covariant field equations, and lead to several additional terms not present at the 
classical level [23]. In other words, the correct field equations for a running G are 
not given simply by a naive replacement G — )■ G(f), which would lead to incorrect 
results, and violate general covariance. 

It is common practice at this point to write an equation for the density contrast 
5{a) as a function not of t, but of the scale factor a{t), by utilizing the identities 

m=aH'-M (75) 

where / is any function of t, and H = a{t)/a{t) is the Hubble constant. This last 
quantity can then be obtained from the zero-th order tt field equation 

3 =8nGop + X (77) 
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re-written in terms of H(a) as 

2 / . X 2 



l + z 



n {i + zf + nR {i + zf + nx 



(78) 



with a = where z is a red shift, Hq the Hubble constant evaluated today, and 
fl the matter (baryonic and dark) density, the space curvature contribution 
corresponding to a curvature k term, and fi^the dark energy part. 



A 



with + JIa = 1 



After introducing the parameter as the cosmological constant fraction 



= 0> 



1 - a 



(79) 



(80) 



one then obtains an equation for the density contrast (5(a) in the normal (i.e. non- 
running G) case 



da? 



dlnHia) 



with growing solution 



da 



da{a) 



d6{a) 
da 



AttGo 



1 



a^H{ay 



■p{a)S{a)^0 (81) 



,1 11 
a-2Fi(-,l;-; 



(82) 



where 2^1 is a hypergeometric function. 

To determine the quantum correction to 5{a) originating from G(n) in Eq. (j54p . 
one sets 



S{a) = <5o(a)[l + co/(a)] 



(83) 



and then uses this linear Ansatz to find the form of /(a) to lowest order in cq. The 
correction is unambiguously determined from the field equations with a running 
Netwon constant G(n), but here the running of G (due to the choice of variables) 
naturally takes on the form 



G(a) = Go 



1 



ao 



(84) 



with a scale factor a « ao corresponding to a mode for which fc « tt/^ (thus ao 
is not necessarily identified with the scale factor "today"). Furthermore one has 
J — 3/21/ with v — 1/3, as determined from lattice gravity in four dimensions. 

What then remains to be done is to compute the growth index / = 9 In (5/9 In a, 
and from it the growth index exponent 7 defined through f — W [24]. Ultimately 
one is interested in the value for this quantities in the vicinity of a current matter 
fraction $7 k, 0.25. For a constant (i.e. not scale dependent) Newton's constant 
one has the well know result / = 0.6028 and exponent 7 = 0.5562. An explicit 
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calculation in the presence of a running G and for a matter fraction w 0.25 gives, 
to lowest linear order in the small quantum correction cq [25], 

7 = 0.5562 -CqCa- (85) 

Here Ca is the coefficient of the quantum correction in the expression for G{a), and 
therefore fixed by the underlying lattice gravity calculations, and Cq an explicitly 
calculable numerical constant that comes out of the solution of the full effective 
covariant field equations for 6 (a). 

The perturbed RW metric is well suited for discussing matter perturbations, but 
occasionally one finds it more convenient to use a different metric parametrization, 
such as the one derived from the conformal Newtonian (cN) gauge line element 

dr^ a^{t){{l + 2i})df - {1 -2(j})5ijdx'dx^} (86) 

with Conformal Newtonian potentials '0(x, t) and (j)ix,t). In this gauge, and in the 
absence of a G{0)^ the unperturbed equations are 

dV Stt 2 - 
a J 6 

^(-)=-yGa^(p-4-3,), (87) 

in the absence of spatial curvature (fc = 0). In the presence of a running G these 
again need to be modified, in accordance with Eqs. (|671) . (|66l) and (|68|) . A cosmolog- 
ical constant can be conveniently included in the p and p, with px = X/8ttG — —p\. 
In this gauge scalar perturbations are characterized by Fourier modes ip{k,t) and 
(j){k,t), and the first order Einstein field equations in the absence of G{0) read [26] 

fc2 + 3 - ((/)+- V I = 471602 6T% 
a \ a J 

^ + = An G {p + p) 

a a" \ , fc2 , , , , 47r „ 2 



a V 



where the perfect fluid energy-momentum tensor is given to linear order in the 
perturbations dp = p ~ p and Sp = p ~ p by 

n ^-ip + dp) 

T\ ^{p + p)v, = -T\ 

T) = {p + 5p) 5) + S'^ E\ = (89) 

and one has allowed for an anisotropic shear perturbation S*^ to the perfect fluid 
form . The two quantities 9 and a are commonly defined by 

ip + p)0 = i F ST"j {p + p)a = -(kSj - (90) 
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with S'^ = T'j — S''jT\/3 the traceless component of T^ . For a perfect fluid 6 is the 
divergence of the fluid velocity, 9 = ik^Vj, with — dx^ /dt the small velocity of 
the fluid. The field equations imply, by consistency, the covariant energy momentum 
conservation law 

6 = -{l + w) {e-3(j))-3- w) 5 

a \dp J 

e = --{l-3w)9- 9 + kH - k^a + k^^P (91) 

a 1 + w 1 + w op 

and relate the matter fields S, a and 9 to the metric perturbations (j) and ip. where 
S is the matter density contrast S = 5p/p, and w is the equation of state parameter 
w — p/p. 

In the presence of a G(n) the above equations need to be re-derived and amended 
[25], starting from the covariant field equations of Eq. (|56l) in the cN gauge of 
Eq. (|86p. with zero-th order modified field equations as in Eqs. (1661) and (j67l) . using 
the expansion for G{0) given in Eq. (j70p . but now in terms of the new cN gauge 
potentials (p and ip. One key question is then the nature of the vacuum-polarization 
induced anisotropic shear perturbation correction S*^- appearing in the covariant 
effective field equations analogous to Eqs. (l88l) . but derived with a G{0). In par- 
ticular one would expect the quantum correction to the energy momentum tensor 
appearing on the r.h.s. of Eq. (|56|) to contribute new terms to the last of Eqs. (f88| . 
which could then account for a non-zero stress a, and thus for a small deviation 
from the classical result for a perfect fiuid, (p = ip. 
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